Genetically Predicted Body Mass Index and Mortality in COPD.

RATIONALE: BMI is associated with COPD mortality, but the underlying mechanisms are unclear. The effect of genetic variants aggregated into a polygenic score may elucidate causal mechanisms and predict risk. OBJECTIVES: To examine the associations of genetically predicted BMI with all-cause and cause-specific mortality in COPD. METHODS: We developed a polygenic score for BMI (PGSBMI) and tested for associations of the PGSBMI with all-cause, respiratory, and cardiovascular mortality in participants with COPD from the COPDGene, ECLIPSE, and Framingham Heart studies. We calculated the difference between measured BMI and PGS-predicted BMI (BMIdiff) and categorized participants into groups of discordantly low (BMIdiff < 20th percentile), concordant (BMIdiff between 20th - 80th percentile), and discordantly high (BMIdiff > 80th percentile) BMI. We applied Cox models, examined potential non-linear associations of the PGSBMI and BMIdiff with mortality, and summarized results with meta-analysis. MEASUREMENTS AND MAIN RESULTS: We observed significant non-linear associations of measured BMI and BMIdiff, but not PGSBMI, with all-cause mortality. In meta-analyses, a one standard deviation increase in the PGSBMI was associated with an increased hazard for cardiovascular mortality (HR=1.29, 95% CI=1.12-1.49), but not with respiratory or all-cause mortality. Compared to participants with concordant measured and genetically predicted BMI, those with discordantly low BMI had higher mortality risk for all-cause (HR=1.57, CI=1.41-1.74) and respiratory death (HR=2.01, CI=1.61-2.51). CONCLUSIONS: In people with COPD, higher genetically predicted BMI is associated with higher cardiovascular mortality but not respiratory mortality. Individuals with discordantly low BMI have higher all-cause and respiratory mortality compared to those with concordant BMI.

Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function.

BACKGROUND: Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function. OBJECTIVES: We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions. METHODS: We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings. RESULTS: Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers. CONCLUSION: Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Key variants via the Alzheimer's Disease Sequencing Project whole genome sequence data.

INTRODUCTION: Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci. METHODS: We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases = 2184, N controls = 2383) and targeted analyses in subpopulations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants. RESULTS: Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses. DISCUSSION: This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Genome-wide association study of preserved ratio impaired spirometry (PRISm).

BACKGROUND: Preserved ratio impaired spirometry (PRISm) is defined as a forced expiratory volume in 1 s (FEV1) <80% predicted and FEV1/forced vital capacity ≥0.70. PRISm is associated with respiratory symptoms and comorbidities. Our objective was to discover novel genetic signals for PRISm and see if they provide insight into the pathogenesis of PRISm and associated comorbidities. METHODS: We undertook a genome-wide association study (GWAS) of PRISm in UK Biobank participants (Stage 1), and selected single nucleotide polymorphisms (SNPs) reaching genome-wide significance for replication in 13 cohorts (Stage 2). A combined meta-analysis of Stage 1 and Stage 2 was done to determine top SNPs. We used cross-trait linkage disequilibrium score regression to estimate genome-wide genetic correlation between PRISm and pulmonary and extrapulmonary traits. Phenome-wide association studies of top SNPs were performed. RESULTS: 22 signals reached significance in the joint meta-analysis, including four signals novel for lung function. A strong genome-wide genetic correlation (rg) between PRISm and spirometric COPD (rg=0.62, p<0.001) was observed, and genetic correlation with type 2 diabetes (rg=0.12, p=0.007). Phenome-wide association studies showed that 18 of 22 signals were associated with diabetic traits and seven with blood pressure traits. CONCLUSION: This is the first GWAS to successfully identify SNPs associated with PRISm. Four of the signals, rs7652391 (nearest gene MECOM), rs9431040 (HLX), rs62018863 (TMEM114) and rs185937162 (HLA-B), have not been described in association with lung function before, demonstrating the utility of using different lung function phenotypes in GWAS. Genetic factors associated with PRISm are strongly correlated with risk of both other lung diseases and extrapulmonary comorbidity.

Dynamic trajectories of life-threatening mass effect in patients with large middle cerebral artery stroke.

Background: Life-threatening, space-occupying mass effect due to cerebral edema and/or hemorrhagic transformation is an early complication of patients with middle cerebral artery (MCA) stroke. Little is known about longitudinal trajectories of laboratory and vital signs leading up to radiographic and clinical deterioration related to this mass effect. Methods: We curated a granular retrospective dataset of 635 patients with large middle cerebral artery (MCA) stroke totaling 108,547 data points for repeated measurements of 10 covariates, and 40 time-independent covariates. We assessed longitudinal trajectories of the 10 longitudinal variables during the 72 hours preceding three outcomes representative of life-threatening mass effect: midline shift (MLS) ≥5mm, pineal gland shift (PGS) >4mm, and decompressive hemicraniectomy (DHC). We used a "backward looking" trajectory approach. Patients were aligned according to the time of outcome occurrence and the trajectory of each variable was assessed prior to that outcome by accounting for both cases and non-cases. Results: Of 635 patients, 49% were female, and mean age was 69 years. Thirty five percent of patients had MLS ≥ 5mm, 24.1% had PGS >4mm, and DHC occurred in 10.7%. For the three outcomes of interest, backward-looking trajectories showed mild increases in white blood cell count (10 up to 11 K/UL within 72 hours), temperature (up to half a degree within 24 hours), and sodium (1-3 mEq/L within 24 hours) leading up to outcomes. We also observed a decrease in heart rate (75 - 65 beats per minute) 24 hours prior to DHC. Conclusions: Univariable longitudinal profiling showed that temperature, white blood cell count, and sodium increase prior to radiographic and clinical indicators of space-occupying mass effect. These findings will inform development of multivariable dynamic risk models to aid prediction of life-threatening space-occupying mass effect.

Bik promotes proteasomal degradation to control low-grade inflammation.

Although chronic low-grade inflammation does not cause immediate clinical symptoms, over the longer term, it can enhance other insults or age-dependent damage to organ systems and thereby contribute to age-related disorders, such as respiratory disorders, heart disease, metabolic disorders, autoimmunity, and cancer. However, the molecular mechanisms governing low-level inflammation are largely unknown. We discovered that Bcl-2-interacting killer (Bik) deficiency causes low-level inflammation even at baseline and the development of spontaneous emphysema in female but not male mice. Similarly, a single nucleotide polymorphism that reduced Bik levels was associated with increased inflammation and enhanced decline in lung function in humans. Transgenic expression of Bik in the airways of Bik-deficient mice inhibited allergen- or LPS-induced lung inflammation and reversed emphysema in female mice. Bik deficiency increased nuclear but not cytosolic p65 levels because Bik, by modifying the BH4 domain of Bcl-2, interacted with regulatory particle non-ATPase 1 (RPN1) and RPN2 and enhanced proteasomal degradation of nuclear proteins. Bik deficiency increased inflammation primarily in females because Bcl-2 and Bik levels were reduced in lung tissues and airway cells of female compared with male mice. Therefore, controlling low-grade inflammation by modifying the unappreciated role of Bik and Bcl-2 in facilitating proteasomal degradation of nuclear proteins may be crucial in treating chronic age-related diseases.

Follow-up ASPECTS improves prediction of potentially lethal malignant edema in patients with large middle cerebral artery stroke.

BACKGROUND: Recent studies have shown that follow-up head CT is a strong predictor of functional outcomes in patients with middle cerebral artery stroke and mechanical thrombectomy. We sought to determine whether total and/or regional follow-up Alberta Stroke Program Early CT Score (ASPECTSfu) are associated with important clinical outcomes during hospitalization and improve the performance of clinical prediction models of potentially lethal malignant edema (PLME). METHODS: We conducted a retrospective study of patients at three medical centers in a major North American metropolitan area with baseline and follow-up head CTs after large middle cerebral artery stroke between 2006 and 2022. We used multivariable logistic regression to test the association of total and regional ASPECTSfu with PLME (cerebral edema related death or surgery), adjusting for total baseline ASPECTS, age, sex, admission glucose, tissue plasminogen activator, and mechanical thrombectomy. We compared existing clinical risk models with and without total or regional ASPECTSfu using area under the curve. RESULTS: In our 560 patient cohort, lower total ASPECTSfu was significantly associated with higher odds of PLME when adjusting for confounders (OR 1.69, 95% CI 1.49 to 2.0), and improved model discrimination compared with existing models and models using baseline ASPECTS. Deep territory involvement (OR 2.46, 95% CI 1.53 to 4.01) and anterior territory involvement (OR 3.23, 95% CI 1.88 to 5.71) were significantly associated with PLME. CONCLUSIONS: Lower ASPECTSfu and certain locations on regional ASPECTSfu, including deep and anterior areas, were significantly associated with PLME. Including ASPECTSfu information improved discrimination of established edema prediction models and could be used immediately to help facilitate clinical management decisions and prognostication.

Frequency of Variants in Mendelian Alzheimer's Disease Genes within the Alzheimer's Disease Sequencing Project (ADSP).

Objective: To characterize the consequences of loss-of-function (LoF), predicted damaging missense (DM), and previously-reported clinically-relevant variants in three Mendelian Alzheimer's disease (AD) genes - presenilin-2 ( PSEN2 ), presenilin-1 ( PSEN1 ), and amyloid precursor protein ( APP ) - within the participants from the Alzheimer's Disease Sequencing Project (ADSP) whole genome sequence (WGS) and whole exome sequence (WES) data. Methods: We identified rare variants (MAF <1%) previously-reported in PSEN2 , PSEN1, and APP in the available ADSP sample of 14,641 individuals with WGS and 16,849 individuals with WES available for research-use (N total = 31,490). We additionally curated variants in these three genes from ClinVar, OMIM, and Alzforum and report carriers of variants in clinical databases as well as LoF and predicted DM variants in these genes. Results: We detected 31 previously-reported clinically-relevant variants with alternate alleles observed within the ADSP: 4 variants in PSEN2 , 25 in PSEN1 , and 2 in APP . Thirty-eight clinical variants with conflicting pathogenicity interpretation within ClinVar or across the databases were identified along with 12 additional LoF and 197 additional DM variants. The overall variant carrier rate for the 31 clinically-relevant variants in the ADSP was 0.3%. We observed 78 individuals carrying at least one clinically-relevant variant, 79.5% were cases compared to 3.9% controls. In those with AD, we observed that the mean age of onset of AD among carriers of these clinically-relevant variants was 19.6 ± 1.4 years earlier compared with non-carriers (p-value=7.8×10 -57 ), and the average age of onset of AD is 5 years earlier in carriers of an additional LoF variant (n=5) compared with non-carriers. Conclusion: The ADSP data permit further characterization of previously-reported AD clinically-relevant variants. A small proportion of individuals in the ADSP are carriers of a previously-reported clinically-relevant variant allele for AD and these participants have significantly earlier age of AD onset compared to non-carriers. Furthermore, we observed additional LoF variants that potentially contribute to clinical presentation of AD.

Multi-tissue epigenetic analysis identifies distinct associations underlying insulin resistance and Alzheimer's disease at CPT1A locus.

BACKGROUND: Insulin resistance (IR) is a major risk factor for Alzheimer's disease (AD) dementia. The mechanisms by which IR predisposes to AD are not well-understood. Epigenetic studies may help identify molecular signatures of IR associated with AD, thus improving our understanding of the biological and regulatory mechanisms linking IR and AD. METHODS: We conducted an epigenome-wide association study of IR, quantified using the homeostatic model assessment of IR (HOMA-IR) and adjusted for body mass index, in 3,167 participants from the Framingham Heart Study (FHS) without type 2 diabetes at the time of blood draw used for methylation measurement. We identified DNA methylation markers associated with IR at the genome-wide level accounting for multiple testing (P < 1.1 × 10-7) and evaluated their association with neurological traits in participants from the FHS (N = 3040) and the Religious Orders Study/Memory and Aging Project (ROSMAP, N = 707). DNA methylation profiles were measured in blood (FHS) or dorsolateral prefrontal cortex (ROSMAP) using the Illumina HumanMethylation450 BeadChip. Linear regressions (ROSMAP) or mixed-effects models accounting for familial relatedness (FHS) adjusted for age, sex, cohort, self-reported race, batch, and cell type proportions were used to assess associations between DNA methylation and neurological traits accounting for multiple testing. RESULTS: We confirmed the strong association of blood DNA methylation with IR at three loci (cg17901584-DHCR24, cg17058475-CPT1A, cg00574958-CPT1A, and cg06500161-ABCG1). In FHS, higher levels of blood DNA methylation at cg00574958 and cg17058475 were both associated with lower IR (P = 2.4 × 10-11 and P = 9.0 × 10-8), larger total brain volumes (P = 0.03 and P = 9.7 × 10-4), and smaller log lateral ventricular volumes (P = 0.07 and P = 0.03). In ROSMAP, higher levels of brain DNA methylation at the same two CPT1A markers were associated with greater risk of cognitive impairment (P = 0.005 and P = 0.02) and higher AD-related indices (CERAD score: P = 5 × 10-4 and 0.001; Braak stage: P = 0.004 and P = 0.01). CONCLUSIONS: Our results suggest potentially distinct epigenetic regulatory mechanisms between peripheral blood and dorsolateral prefrontal cortex tissues underlying IR and AD at CPT1A locus.

Heterogeneous effects on type 2 diabetes and cardiovascular outcomes of genetic variants and traits associated with fasting insulin.

Hyperinsulinemia is a complex and heterogeneous phenotype that characterizes molecular alterations that precede the development of type 2 diabetes (T2D). It results from a complex combination of molecular processes, including insulin secretion and insulin sensitivity, that differ between individuals. To better understand the physiology of hyperinsulinemia and ultimately T2D, we implemented a genetic approach grouping fasting insulin (FI)-associated genetic variants based on their molecular and phenotypic similarities. We identified seven distinctive genetic clusters representing different physiologic mechanisms leading to rising FI levels, ranging from clusters of variants with effects on increased FI, but without increased risk of T2D (non-diabetogenic hyperinsulinemia), to clusters of variants that increase FI and T2D risk with demonstrated strong effects on body fat distribution, liver, lipid, and inflammatory processes (diabetogenic hyperinsulinemia). We generated cluster-specific polygenic scores in 1,104,258 individuals from five multi-ancestry cohorts to show that the clusters differed in associations with cardiometabolic traits. Among clusters characterized by non-diabetogenic hyperinsulinemia, there was both increased and decreased risk of coronary artery disease despite the non-increased risk of T2D. Similarly, the clusters characterized by diabetogenic hyperinsulinemia were associated with an increased risk of T2D, yet had differing risks of cardiovascular conditions, including coronary artery disease, myocardial infarction, and stroke. The strongest cluster-T2D associations were observed with the same direction of effect in non-Hispanic Black, Hispanic, non-Hispanic White, and non-Hispanic East Asian populations. These genetic clusters provide important insights into granular metabolic processes underlying the physiology of hyperinsulinemia, notably highlighting specific processes that decouple increasing FI levels from T2D and cardiovascular risk. Our findings suggest that increasing FI levels are not invariably associated with adverse cardiometabolic outcomes.

Association Between Whole Blood-Derived Mitochondrial DNA Copy Number, Low-Density Lipoprotein Cholesterol, and Cardiovascular Disease Risk.

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross-sectional and prospective association analyses of blood-derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole-genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P<0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1-SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04-1.12; P<0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (ß=0.091; P=0.11) or in the reverse direction (ß=-0.012; P=0.076). Additional bidirectional Mendelian randomization analyses revealed that low-density lipoprotein cholesterol had a causal effect on mtDNA CN (ß=-0.084; P<0.001), but the reverse direction was not significant (P=0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low-density lipoprotein cholesterol level (P=0.52), whereas there was a strong direct causal effect of higher low-density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (ß=-0.092; P<0.001). Conclusions Our findings indicate that high low-density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

A Type 1 Diabetes Polygenic Score Is Not Associated With Prevalent Type 2 Diabetes in Large Population Studies.

Context: Both type 1 diabetes (T1D) and type 2 diabetes (T2D) have significant genetic contributions to risk and understanding their overlap can offer clinical insight. Objective: We examined whether a T1D polygenic score (PS) was associated with a diagnosis of T2D in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium. Methods: We constructed a T1D PS using 79 known single nucleotide polymorphisms associated with T1D risk. We analyzed 13 792 T2D cases and 14 169 controls from CHARGE cohorts to determine the association between the T1D PS and T2D prevalence. We validated findings in an independent sample of 2256 T2D cases and 27 052 controls from the Mass General Brigham Biobank (MGB Biobank). As secondary analyses in 5228 T2D cases from CHARGE, we used multivariable regression models to assess the association of the T1D PS with clinical outcomes associated with T1D. Results: The T1D PS was not associated with T2D both in CHARGE (P = .15) and in the MGB Biobank (P = .87). The partitioned human leukocyte antigens only PS was associated with T2D in CHARGE (OR 1.02 per 1 SD increase in PS, 95% CI 1.01-1.03, P = .006) but not in the MGB Biobank. The T1D PS was weakly associated with insulin use (OR 1.007, 95% CI 1.001-1.012, P = .03) in CHARGE T2D cases but not with other outcomes. Conclusion: In large biobank samples, a common variant PS for T1D was not consistently associated with prevalent T2D. However, possible heterogeneity in T2D cannot be ruled out and future studies are needed do subphenotyping.

Bronchial epithelial gene expression and interstitial lung abnormalities.

INTRODUCTION: Interstitial lung abnormalities (ILA) often represent early fibrotic changes that can portend a progressive fibrotic phenotype. In particular, the fibrotic subtype of ILA is associated with increased mortality and rapid decline in lung function. Understanding the differential gene expression that occurs in the lungs of participants with fibrotic ILA may provide insight into development of a useful biomarker for early detection and therapeutic targets for progressive pulmonary fibrosis. METHODS: Measures of ILA and gene expression data were available in 213 participants in the Detection of Early Lung Cancer Among Military Personnel (DECAMP1 and DECAMP2) cohorts. ILA was defined using Fleischner Society guidelines and determined by sequential reading of computed tomography (CT) scans. Primary analysis focused on comparing gene expression in ILA with usual interstitial pneumonia (UIP) pattern with those with no ILA. RESULTS: ILA was present in 51 (24%) participants, of which 16 (7%) were subtyped as ILA with a UIP pattern. One gene, pro platelet basic protein (PPBP) and seventeen pathways (e.g. TNF-α signalling) were significantly differentially expressed between those with a probable or definite UIP pattern of ILA compared to those without ILA. 16 of these 17 pathways, but no individual gene, met significance when comparing those with ILA to those without ILA. CONCLUSION: Our study demonstrates that abnormal inflammatory processes are apparent in the bronchial airway gene expression profiles of smokers with and without lung cancer with ILA. Future studies with larger and more diverse populations will be needed to confirm these findings.

Key variants via Alzheimer's Disease Sequencing Project whole genome sequence data.

INTRODUCTION: Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci. METHODS: We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases=2,184, N controls=2,383) and targeted analyses in sub-populations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants. RESULTS: Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses. DISCUSSION: This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Association of Common and Rare Variants with Alzheimer's Disease in over 13,000 Diverse Individuals with Whole-Genome Sequencing from the Alzheimer's Disease Sequencing Project.

Alzheimer's Disease (AD) is a common disorder of the elderly that is both highly heritable and genetically heterogeneous. Here, we investigated the association between AD and both common variants and aggregates of rare coding and noncoding variants in 13,371 individuals of diverse ancestry with whole genome sequence (WGS) data. Pooled-population analyses identified genetic variants in or near APOE, BIN1, and LINC00320 significantly associated with AD (p < 5×10-8). Population-specific analyses identified a haplotype on chromosome 14 including PSEN1 associated with AD in Hispanics, further supported by aggregate testing of rare coding and noncoding variants in this region. Finally, we observed suggestive associations (p < 5×10-5) of aggregates of rare coding rare variants in ABCA7 among non-Hispanic Whites (p=5.4×10-6), and rare noncoding variants in the promoter of TOMM40 distinct of APOE in pooled-population analyses (p=7.2×10-8). Complementary pooled-population and population-specific analyses offered unique insights into the genetic architecture of AD.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Clonal Hematopoiesis of Indeterminate Potential (CHIP) and Incident Type 2 Diabetes Risk.

OBJECTIVE: Clonal hematopoiesis of indeterminate potential (CHIP) is an aging-related accumulation of somatic mutations in hematopoietic stem cells, leading to clonal expansion. CHIP presence has been implicated in atherosclerotic coronary heart disease (CHD) and all-cause mortality, but its association with incident type 2 diabetes (T2D) is unknown. We hypothesized that CHIP is associated with elevated risk of T2D. RESEARCH DESIGN AND METHODS: CHIP was derived from whole-genome sequencing of blood DNA in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) prospective cohorts. We performed analysis for 17,637 participants from six cohorts, without prior T2D, cardiovascular disease, or cancer. We evaluated baseline CHIP versus no CHIP prevalence with incident T2D, including associations with DNMT3A, TET2, ASXL1, JAK2, and TP53 variants. We estimated multivariable-adjusted hazard ratios (HRs) and 95% CIs with adjustment for age, sex, BMI, smoking, alcohol, education, self-reported race/ethnicity, and combined cohorts' estimates via fixed-effects meta-analysis. RESULTS: Mean (SD) age was 63.4 (11.5) years, 76% were female, and CHIP prevalence was 6.0% (n = 1,055) at baseline. T2D was diagnosed in n = 2,467 over mean follow-up of 9.8 years. Participants with CHIP had 23% (CI 1.04, 1.45) higher risk of T2D than those with no CHIP. Specifically, higher risk was for TET2 (HR 1.48; CI 1.05, 2.08) and ASXL1 (HR 1.76; CI 1.03, 2.99) mutations; DNMT3A was nonsignificant (HR 1.15; CI 0.93, 1.43). Statistical power was limited for JAK2 and TP53 analyses. CONCLUSIONS: CHIP was associated with higher incidence of T2D. CHIP mutations located on genes implicated in CHD and mortality were also related to T2D, suggesting shared aging-related pathology.

Investigating Associations of Omega-3 Fatty Acids, Lung Function Decline, and Airway Obstruction.

Rationale: Inflammation contributes to lung function decline and the development of chronic obstructive pulmonary disease. Omega-3 fatty acids have antiinflammatory properties and may benefit lung health. Objectives: To investigate associations of omega-3 fatty acids with lung function decline and incident airway obstruction in a diverse sample of adults from general-population cohorts. Methods: Complementary study designs: 1) longitudinal study of plasma phospholipid omega-3 fatty acids and repeated FEV1 and FVC measures in the NHLBI Pooled Cohorts Study and 2) two-sample Mendelian randomization (MR) study of genetically predicted omega-3 fatty acids and lung function parameters. Measurements and Main Results: The longitudinal study found that higher omega-3 fatty acid levels were associated with attenuated lung function decline in 15,063 participants, with the largest effect sizes for the most metabolically downstream omega-3 fatty acid, docosahexaenoic acid (DHA). An increase in DHA of 1% of total fatty acids was associated with attenuations of 1.4 ml/yr for FEV1 (95% confidence interval [CI], 1.1-1.8) and 2.0 ml/yr for FVC (95% CI, 1.6-2.4) and a 7% lower incidence of spirometry-defined airway obstruction (95% CI, 0.89-0.97). DHA associations persisted across sexes and smoking histories and in Black, White, and Hispanic participants, with associations of the largest magnitude in former smokers and Hispanic participants. The MR study showed similar trends toward positive associations of genetically predicted downstream omega-3 fatty acids with FEV1 and FVC. Conclusions: The longitudinal and MR studies provide evidence supporting beneficial effects of higher levels of downstream omega-3 fatty acids, especially DHA, on lung health.

A Polygenic Risk Score for Idiopathic Pulmonary Fibrosis and Interstitial Lung Abnormalities.

Rationale: In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives: We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods: We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results: We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case-control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10-95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10-87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions: A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Discovering cellular programs of intrinsic and extrinsic drivers of metabolic traits using LipocyteProfiler.

A primary obstacle in translating genetic associations with disease into therapeutic strategies is elucidating the cellular programs affected by genetic risk variants and effector genes. Here, we introduce LipocyteProfiler, a cardiometabolic-disease-oriented high-content image-based profiling tool that enables evaluation of thousands of morphological and cellular profiles that can be systematically linked to genes and genetic variants relevant to cardiometabolic disease. We show that LipocyteProfiler allows surveillance of diverse cellular programs by generating rich context- and process-specific cellular profiles across hepatocyte and adipocyte cell-state transitions. We use LipocyteProfiler to identify known and novel cellular mechanisms altered by polygenic risk of metabolic disease, including insulin resistance, fat distribution, and the polygenic contribution to lipodystrophy. LipocyteProfiler paves the way for large-scale forward and reverse deep phenotypic profiling in lipocytes and provides a framework for the unbiased identification of causal relationships between genetic variants and cellular programs relevant to human disease.